Anomalous pairing vibration in neutron-rich Sn isotopes beyond the N = 82 magic 
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Two-neutron transfer associated with the pair correlation in superfluid neutron-rich nuclei is stud- 
ied with focus on low-lying + states in Sn isotopes beyond the N — 82 magic number. We describe 
microscopically the two-neutron addition and removal transitions by means of the Skyrme-Hartree- 
Fock-Bogoliubov mean-field model and the continuum quasiparticie random phase approximation 
formulated in the coordinate space representation. It is found that the pair transfer strength for 
the transitions between the ground states becomes significantly large for the isotopes with A > 140, 
reflecting very small neutron separation energy and long tails of the weakly bound 3p orbits. In 
132_140 Sn, a peculiar feature of the pair transfer is seen in transitions to low-lying excited + states. 
They can be regarded as a novel kind of pair vibrational mode which is characterized by an anoma- 
lously long tail of the transition density extending to far outside of the nuclear surface, and a large 
strength comparable to that of the ground-state transitions. The presence of the weakly bound 
neutron orbits plays a central role for these anomalous behaviors. 

PACS numbers: 21.10.Pc, 21.10.Re, 21.60.Jz, 25.40.Hs, 27.60.+j 



I. INTRODUCTION 

Recently two-neutron transfer processes have attracted 
renewed interests thanks to the opportunities of two- 
neutron transfer reaction experiments using the beams of 
radioactive neutron- rich isotopes [lH5j . In the past a con- 
siderable number of experimental and theoretical studies 
of two-neutron transfer have been accumulated for sta- 
ble nuclei [6T4l2l]. One of the central concepts established 
there is that the two-neutron transfer amplitude is influ- 
enced by collective excitation modes caused by the su- 
perfluidity or the pair correlation [Hj]. In the case 
of the open-shell nuclei, the strong two-neutron transfer 
transitions connecting the ground states emerge; it is a 
Nambu-Goldstonc mode related to the rotational sym- 
metry with respect to the phase of the pair condensate, 
called the pairing rotation. The pairing vibration, which 
corresponds to a vibration of the pairing gap, is another 
collective mode producing a low-lying excited + state, 
but the intensity of the associated pair transfer is signif- 
icantly smaller than the pairing rotation. 

When neutron-rich nuclei are concerned, the above 
conventional picture of the pairing collectivity may be 
modified. We can expect this because neutron-rich nuclei 
often accompany skin or halo, i.e., low density distribu- 
tions of weakly bound neutrons surrounding the nucleus, 
and also because the neutron pairing in neutron matter 
is predicted to be stron g or close to the strong-coupling 
regime at low densities |13l - [l3 |. The expected enhance- 
ment of the pair correlation may result in unusual proper- 
ties in two-neutron transfers. Dobaczewski et al.[20| has 
pointed out the possibility that the two-neutron transfer 
associated with the pairing rotation may be enhanced in 
neutron-rich nuclei because of the surface enhanced pair- 
ing. Influences of the surface enhancement on the pair 
vibrational two-neutron transfer populating excited + 



states and the first 2 + states in neutron-rich Sn isotopes 
are discussed in Refs. [HI, H2| and in Ref. [23j], respec- 
tively. Recent experiments of (Pjt) reaction on light-mass 
neutron-rich nuclei such as n Li (HHIll and M He [3,[1- 
[BJ point to the importance of the pair correlation in these 
typical halo or skin nuclei. The giant pairin g vi bration 
in neutron-rich nuclei has also been discussed 26, 271. 



In the present paper we study the pairing collective 
modes in heavy-mass neutron-rich nuclei in order to ex- 
plore new features in the pairing rotation, the pairing 
vibration, and the associated two-neutron transfer am- 
plitudes. We focus on excitation modes with monopole 
multipolarity (L = 0), for which the pairing collectivity 
is expected most strongly. As a target of the analysis, 
we choose the Sn isotopic chain ranging from the proton- 
rich side 100 Sn (N = 50) to very neutron-rich isotopes 
with j4 ^ 150 beyond the N — 82 magic number ( 132 Sn 
). Our theoretical tool is the Hartree-Fock-Bogoliubov 
mean-field model using the Skyrme energy density func- 
tional and the continuum quasiparticie random phase ap- 
proximation, and it is the same as adopted in Ref. [231] . 
A similar approach is made in Khan et al. [2l|, [26| , where 
the L = pairing vibration in neutron-rich Sn isotopes is 
also studied. Compared with Ref.[2l|, we perform more 
systematic and detailed analyses, which eventually bring 
us a new finding and some differences. Among all, the 
present analysis predicts a novel type of the pairing vi- 
bration that emerges in 134_140 Sn. 
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II. SKYRME-HARTREE-FOCK-BOGOLIUBOV 
MEAN-FIELD PLUS QRPA APPROACH 

A. The model 

We describe the neutron pair correlation by means 
of the Hartree-Fock-Bogoliubov (HFB) mean-field 
theory [28( and the continuum quasiparticle random 
phase approximation (QRPA) [2^, We assume the 

spherical symmetry of the mean-fields and densities as- 
sociated with the ground states of the subshell-closed Sn 
isotopes. The model adopted in the present study is the 
same as that in Ref . [23j . 

The starting point of the model is the energy density 
functional, which we construct from the Skyrme interac- 
tion with the parameter set SLy4[3l|, and the density- 
dependent delta interaction (DDDI) adopted as an effec- 
tive pairing force [H GJ, IM lH-[4l|. The DDDI is given 
by 

vf(r,r') = ^(l-P a )V q (r)S(r-r') (q = n,p) (1) 

where V q (r) is the pairing interaction strength and is a 
function of the neutron and proton densities. We adopt 
the following form 



V n (r) =v 1 - T) 



with the parameters vo 



Pn(r) 



-458.4 MeV fm 3 , p c 



(2) 



0.08fnT 



, a 



0.59, and r) = 0.71 [16|, |2jJ |39| . Here the 



parameter vq representing the strength in the free space 
is chosen to reproduce the scattering length a = —18.5 
fm of the bare nn-interaction in the 1 S channel, and a is 
to fit the density-dependent BCS pairing gap of neutron 
matter. The parameter rj is adjusted to reproduce the 
experimental pair gap in 120 Sn. We call this parameter 
set DDDI-bare' as it is motivated by the bare nuclear 
force. This interaction has a strong density dependence. 
For comparison, we also use the volume pairing interac- 
tion (vq = —195 MeV fm 3 and ij = 0) with no density- 
dependence, and also the mix pairing interaction[36, 37] 
( Vq = —292 MeV fm 3 ) with intermediate density depen- 
dence. See Ref. |23[ for details. The parameter vq of the 
volume and the mix pairing interactions is determined to 
reproduce the pair gap in 120 Sn. 

The HFB equation is solved in the coordinate space 
representation using the polar coordinate system. The 
radial coordinate space is truncated at r max = 20 fm. 
The density and the pair density of neutrons or protons 
are given as sums of the contributions from the quasipar- 
ticle states, which we truncate with respect to the angu- 
lar momentum quantum numbers and the quasiparticle 
energy. We actually set the maximum orbital angular 
quantum number l max — 12 and the maximum quasipar- 
ticle energy E max — 60 MeV. 

The excited states are described by means of the QRPA 
formulated in the coordinate space representation and 



on the basis of the self-consistent HFB solution. The 
residual interactions to be used in the QRPA are de- 
rived from the Skyrme energy density functional and 
the DDDI, but the Landau-Migdal approximation is em- 
ployed for the residual interaction in the particle-hole 
channel. We impose the outgoing-wave boundary condi- 
tion on the continuum quasiparticle states as described in 
Ref. [29l l3(i| in the case of the isotopes with A > 132 while 
the box boundary condition is adopted [231] for A < 132 
where the neutron Fermi energy is deeper than —6.8 
MeV. In the QRPA response functions, the smearing pa- 
rameter e = 50 keV is introduced so that the obtained 
strength function is convoluted with the Lorentzian with 
the FWHM= 2e = 100 keV. 



B. Strength function and transition density for 
pair transfer modes 

In the present work, we define the strength of the two- 
neutron transfer using matrix elements of the creation 
and annihilation operators of a S = pair of neutrons 
with the angular momentum L: 



LM 



drY LM (f)f(r)^(ri)^(rt), (3) 



Plm = j drY£ M (f)f(r)%(r^ q (ri). (4) 

When we consider transitions from the ground state of a 
nucleus with even N to states in the neighboring nucleus 
with N + 2 through the addition of a neutron pair, we 
evaluate the strength function P^ M 

S PadL (E) =Y,*{E - E iL )\{^ iLM \Pt M \^ )\\ (5) 

iM 

and for transitions via the removal of a neutron pair we 
calculate the strength function 

S PrmL (E) = J2KE- E iL )\(^ iL ^m\Plai\^o)\ 2 . (6) 

iM 

Here \l/o is the ground state with N while %jj is a state 
with the angular quantum numbers LM and the neutron 
number N±2, and En, is the transition energy measured 
from the ground state of the residue nucleus with N ± 2. 

For a specific transition to a given final state popu- 
lated via the pair-addition and removal operators, the 
transition densities 

Pt d) (r) = (y iL M\^(rl)^(rt)\*o) 

= Yl M (r)Pt\r), (7) 

if m) (r) = (¥iLwWr W(r40|*o) 

= Yl M (r)P^ m \r) (8) 

may be defined. We also evaluate the reduced transition 
probabilities defined by 



B(PadL; gs - 



i) = £|<vW|p2mI*o}| 2 

M 



3 



= (2L + 1) 



r 2 Pt\r)dr 



B(PrmL;gs^i) = ^ K^mI^ 



(2L + 1) 



r 2p(rm) 



(r)dr 



(9) 



(10) 



For the procedure to evaluate the strength functions, the 
transition densities and the transition strengths in the 
QRPA scheme, we refer the readers to Refs. [23l. l29l. [30| . 



III. GROUND STATE PROPERTIES AND 
PAIRING ROTATION 

Let us first discuss the pairing rotation, i.e., the 
monopole transition with L = from the ground state 
of a supcrfluid open-shell even- TV isotopes to the ground 
state of neighboring N ± 2 isotopes. The transition den- 
sity and the transition matrix element can be calculated 
at the level of the static HFB as the transition density 
for this transition can be approximated as 



where p~(r) is the pair density defined by 

p(r) = (y \Yyi(ra)^{ra)\* ) 



\p{r) 



(11) 



(12) 



The radial transition density is given as 

PS&(T) = v^P~(r). (13) 

In this approximation, there is no distinction between 
the addition and removal modes. The strength of the 
ground-state transfer is calculated as 

2 



£?(Pad/rmO; gs — > gs) = 



ir r p(r)dr 



(14) 



Note that in the closed-shell isotopes 100 Sn and 132 Sn 
we calculate the strength and the transition density in 
the QRPA instead of Eqs. flU) and ((14|) which are not 
appropriate to nuclei with vanishing pairing gap. 

Figure HJa) is the plot of the calculated strength 
_B(Pad/rmO) of the ground-state two-neutron transfer. 
The strength is enhanced significantly as the mass num- 
ber (the neutron number) exceeds A = 140 (N = 90). 
For the isotopes with A > 140, the absolute magnitude 
of the pair transfer strength reaches more than twice the 
maximum strength in the region 100 < A < 132. In 
the conventional BCS approximation, the pair transfer 
strength is proportional to (A/G) 2 oc A 2 where A is the 
pairing gap and G is the force strength of the seniority 
pairing force @, H, 0|. In Fig. Hfb), we plot the square 
A 2 ^ of the average pairing gap 



A n7 ; 



J p(r)A(r)dr 
J p(r)dr 



(15) 



of neutrons. Comparing the isotopic trends of the 
two quantities, we see relative enhancement of the two- 
neutron transfer strength by a factor of two in A > 140, 
and also in 132 < A < 140. For the latter isotopes, the 
absolute magnitude of the average pairing gap is small, 
but the two-neutron transfer strength is comparable to 
that of more stable isotopes in the region 100 < A < 132. 

The origin of the enhancement is clarified by looking 
at the pair transition density, which is shown in Fig. [Ha) 
and (b) for 120 < A < 132 and 132 < A < 150, re- 
spectively. It is seen that the profile of the transition 
density suddenly changes as the neutron number exceeds 
the N = 82 magic number and N = 90. The transition 
density for 132 < A < 150 extends outside the surface, 
reaching r ~ 11 fm for 132 < A < 140, and r ~ 14 fm 
for 140 < A < 150. The amplitude in the exterior re- 
gion r ^ 7 fm for A > 132 is evidently larger than those 
for A < 132, where the amplitudes extend only up to 
r ~ 9 fm. Comparing the results for A — 120 and for 
A = 144, for instance, the maximum values of the am- 
plitude around the nuclear surface r ~ 6 fm are approxi- 
mately the same, but because of the large spatial exten- 
sion of the transition density, the pair transfer strength 
in 144 Sn is larger by a factor of ~ 2 (cf. FigfTJ). 

The reason for the spatial extension of the pair tran- 
sition density to develop suddenly beyond TV = 82 and 
N = 90 can be ascribed to the shell gap at N = 82 
and properties of the neutron single-particle states. We 
here note that the transition density of the pair rota- 
tional mode, i.e., the pair density p(r) is written as 
a coherent sum of contributions of quasiparticle states, 
and the quasiparticle states with lower excitation energy 
(i.e., those originating from orbits close to the Fermi en- 
ergy) have larger contributions. The calculated Hartree- 
Fock single-particle energies for neutrons in 132 Sn are 
e HF = -1.99,-0.25,0.26 MeV for the 2/ 7/2 , 3p 3/2 , and 
3pi/2 orbits located above the N — 82 gap, respectively ( 
3pi/2 is an unbound resonance), and the hxi/2 orbit with 
ghf = —7.68 MeV is located below the shell gap. For 
the 132 < A < 140 isotopes (where the neutron Fermi 
energy is located near the position of 2/ 7 / 2 ), the main 
component of the transition density originates from this 
orbit. Since the binding energy of 2/ 7 / 2 is rather small, 
the tail of its wave function extends to outside, leading 
to the long tail in the pair transition density. When the 
neutron number exceeds N = 90 (A = 140), the next sin- 
gle particle orbits 3p3/ 2 and 3pi/ 2 give large contribution 
to the pair density. Since these p orbits have very small 
binding or are unbound, the spatial extension further de- 
velops in isotopes with N > 90 (A > 140). 



IV. PAIRING VIBRATION 

A. Strength function 

We now discuss the two-neutron transfer modes popu- 
lating excited + states. Figure |3J a), (b),(c) and (d) show 
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FIG. 1: (a) The neutron pair transfer strength B(Pad/rmO; gs — ¥ gs) of the ground-state transitions calculated for the even-even 
Sn isotopes. The diamonds connected with the solid line are the results obtained with the pairing interaction DDDI-bare' while 
the circles and the crosses are those with the mix and the volume pairing interactions, respectively. (b)The squared average 
pairing gap A 2 „ of neutrons. 
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FIG. 2: (a) The neutron pair transition density r 2 P^^J^\r) for the ground-state transitions in the isotopes 120 130 Sn and 
134 Sn. (b) The same as (a), but in 132_140 Sn and 142_150 Sn. The transition density for 132 Sn is the pair-addition transition 
density r 2 P^^i (r) calculated in the QRPA. 



the strength functions Sp a do{E) and Sp rm o(E) for the 
monopole pair-addition and -removal modes in four rep- 
resentative isotopes 102 Sn, 120 Sn, 134 Sn and 142 Sn. Here, 
120 Sn is a mid-shell isotope near the stability line. 134 Sn 
is a representative of neutron-rich unstable isotope which 
is located just beyond the neutron magic number N = 82 
while 102 Sn has the same closed-shell plus two-particle 
configuration, but at the neutron-deficient side. In the 
most neutron-rich isotope 142 Sn, the neutron Fermi en- 
ergy is located around the weakly bound 3p3/2 orbit. 

All the four isotopes exhibit low-lying peaks in the ex- 
citation energy range of 2 E 4 MeV, and these peaks 
may be regarded as the pairing vibrational modes. They 
corresponds to the lowest QRPA solutions in the cases 
of 102 Sn, 120 Sn and 134 Sn, but second lowest solution in 
the case of 142 Sn. The excitation energies of the pair 
vibrational states are 3.55, 2.40, 3.73 and 2.35 MeV in 
102 Sn, 120 Sn, 134 Sn and 142 Sn, respectively. They are not 
very different from the typical energy E 2A (~ 2MeV) 
expected in the standard pair vibrational model@, HKHj]. 



The most prominent feature seen in Fig(3] is that the 
pair-addition strength associated with the low-lying pair 
vibrational mode in 134 Sn is several times larger than 
those in the other cases. The pair-addition strength of 
the pairing vibration is -B(PadO) = 3.16 in 134 Sn while 
it is 1.18 and 0.40 in 102 Sn and 120 Sn. The magnitude 
of the pair addition strength in 134 Sn is so large that 
it is comparable with the strength £?(Pad/rm0) = 3.61 
associated with the ground-state transition (the pairing 
rotation). The pair- removal strength is negligibly small, 
and hence it is essentially a pure pair-addition mode. We 
also notice that the pair vibrational mode in 134 Sn (at 
E = 3.73 MeV) is located above the one-neutron separa- 
tion energy Si n = 3.31 MeV (indicated with an arrow in 
the figure). We estimate the resonance width being less 
than 1 keV by evaluating the FWHM of the peak mi- 
nus the smoothing width 2e = 100 keV. The very small 
width indicates that the pair vibrational state is a narrow 
resonance even though it is embedded in the continuum. 
The above characteristics indicate clearly that the pair 
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FIG. 3: The strength function Sp a do(E) for the pair-addition mode (plotted in the upper panel with the solid curve) and the 
strength function £p r mo(-E) for the pair-removal mode (the dashed curve in the lower panel) in (a) 102 Sn, (b) 120 Sn, (c) 134 Sn 
and (d) 142 Sn. The arrows indicate the calculated one- and two- neutron separation energies Sin and S'2n- 



vibrational mode in Sn deviates from the conventional 
picture of the pairing vibration. 

The pair vibrational mode in 142 Sn has a character dif- 
ferent from that in 134 Sn as it has the large pair-removal 
strength instead of the pair-addition strength. 

Looking at the strength in 134 Sn and 142 Sn at higher 
excitation energies, it is seen that there exists a smooth 
distribution of the pair-addition strength above the two- 
neutron separation energy 6*2™ ■ At any excitation en- 
ergy above S^n, it is always possible to put two neutrons 
outside the nucleus although these two neutrons imme- 
diately escape out of the initial position. The smooth 
distribution can be ascribed to this process. Since it does 
not correspond to " transfer" reaction, we do not analyze 
below. Concerning the pair-removal strength function, 
sharp and large peaks located around E = 9.1 and 8.6 
MeV are observed in Fig.[3ja) and (c), respectively. This 
is the so called giant pairing vibration [TTi |43T-|451] : it is a 
collective mode of removing two neutrons which occupy 
the single- particle orbits in the next major-shell below 
the valence shell, i.e. the orbits in the N = 50 — 82 shell 
in the case of 134 Sn. In the cases of 120 Sn and 142 Sn, 
similar peaks exist, but at slightly higher energy, and are 
not seen in Fig- EJb) and (d). 



B. Transition density 

We show in Fig|4]the pair-addition and -removal tran- 
sition densities Pi^ {r) and f/™o( r )' an< ^ tne P ar ticle- 
hole transition density p^ =0 (r) for neutrons, associated 
with the low-lying pair vibrational modes in (a) 102 Sn, 
(b) 120 Sn, (c) 134 Sn and (d) 142 Sn. The volume element 
r 2 is multiplied to these quantities in the plots. 

A peculiar feature of the low-lying pair vibrational 
mode in 134 Sn is clearly visible in Fig|4]when it is com- 
pared with the results for 102 Sn and 120 Sn. In the cases of 
102 Sn and 120 Sn, all of the three kinds of transition den- 
sity have increased amplitudes around the nuclear sur- 
face, and the amplitudes diminish quickly far outside the 
surface, r^R lms + 3 fm (r^8 fm). In 134 Sn, however, 

the pair-addition transition density P^=o ( r ) nas a signif- 
icant amplitude even at r = i? rms + 3 fm (r»8 fm) and a 
long tail extends up to r w 15 fm. It is interesting to com- 
pare 102 Sn and 134 Sn (FigQJa) vs. (c)) since both cases 
commonly have the closed-shell plus two-neutron config- 
uration. Although the maximum amplitudes in the two 
cases are about the same, the transition density in 134 Sn 
significantly extends toward outside, resulting in the no- 
ticeable increase of the pair- addition transfer strength in 
134 Sn, i.e., -B(PadO) = 3.16, which is larger by a factor of 
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FIG. 4: (a) The neutron transition densities associated with the pair vibrational mode at E = 3.55 MeV in 102 Sn. The solid, 
the dashed and the dotted curves represent the pair-addition transition density r 2 P i ^=o( r )i the pair-removal transition density 



" 2 f 4 £™o( r )' an< l the particle-hole transition density r 2 p?^_ (r), respectively. (b)(c)(d) The same as (a) but for 120 Sn, 134 Sn 



and Sn, respectively. The position of the arrow indicate the the root mean square radius R r 
density, which are 4.41, 4.67, 4.84 and 4.98 fm in 102 Sn, 120 Sn, 134 Sn and 142 Sn, respectively. 



(r 2 ) of the total nucleon 



three than B(Pad0) = 1.18 in 102 Sn. We thus find that 
the spatial extension of the transition amplitude plays a 
central role to enhance the pair-addition transfer strength 
B(Pad0) in 134 Sn. 

The low- lying pair vibrational mode in 142 Sn exhibits 
also an extended profile in the transition density, but 
in this case the dominant amplitude is the pair-removal 
transition density ^£=o( r ) an< ^ ^ s spatial extension is 
smaller (observed up to r ~ 11 fm) than that of Pjj^n (r) 
in 134 Sn. 



C. Microscopic origin 

In order to clarify the microscopic origin of the anoma- 
lous pair vibrational mode in 134 Sn, we show in Fig. [5](a) 
effects of the residual pairing interaction on the pair- 
addition strength function Sp a do(E). Plotted here is the 
strength function associated with unperturbed neutron 
two-quasiparticle excitations, i.e., the strength obtained 
by neglecting the residual interaction, and it is compared 
with the result of the full QRPA calculation. The first 
peak at E « 1.5MeV and the second one at E w 4.5 MeV 



correspond to the two-quasincutron excitations [2/V/2] 2 
and [3p 3 / 2 ] 2 , respectively. The strengths of these peaks 
are B(Pad0) = 0.18 and 0.15 for [2/ 7/2 ] 2 and [3p 3/2 ] 2 , 
respectively, and they are comparable to a single-particle 
estimate B s . p .(Pad0) = (2j + l)/87t. The large low-lying 
strength does not show up in the unperturbed strength, 
indicating that the collective configuration mixing caused 
by the residual interaction plays a crucial role to bring 
about the large strength of the pair vibrational mode in 
the full QRPA calculation. The collective enhancement 
is a factor of twenty as is estimated from the ratio of the 
strengths (3.16 vs. 0.18, 0.15). 

Figure EJb) shows the pair-addition transition densi- 
ties Pjx^r,{r) associated with unperturbed neutron two- 
quasiparticle excitations [2/V/2] 2 and [3p 3 /2] 2 . It is seen 
that the transition density of the [3p 3 /2] 2 configuration 
has a long tail, extending up to r ~ 15 fm. We here note 
that the quasiparticle energy of the neutron 2/ 7 / 2 and 
3^3/2 states in 134 Sn are 0.75 and 2.25 MeV, and the en- 
ergies of the corresponding Hartree-Fock single-particle 
orbits are enF = —2.14 and —0.36 MeV, respectively. 
Here the binding energy of the 3p 3 / 2 orbit is only one 
fifth of that of 2/7/2, and is smaller than those of the 
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FIG. 5: (a) The pair-addition strength function Sp a do{E) in 134 Sn, calculated for the unperturbed neutron two-quasiparticle 



excitations (plotted with the sold curve) and in the full QRPA (dotted curve), (b) The transition density r 2 P^=o( r ) °f the 



unperturbed neutron two-quasiparticle excitations [2/ 7 / 2 ] (solid curve) and [3p 3 / 2 ] (dashed curve) m 
the transition density of the low- lying pair vibrational mode obtained in the full QRPA (dotted curve). 
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FIG. 6: (a) The pair-removal strength function Sp rm o(i?) in 142 Sn, calculated for the unperturbed neutron two-quasiparticle 
excitations (solid curve) and in the full QRPA (dotted curve), (b) The transition density r 2 P^^ (r) of the unperturbed neutron 



two-quasiparticle excitations [2/7/2] (solid curve) and [3P3/2] (dashed curve) in 
of the low-lying pair vibrational mode obtained in the full QRPA (dotted curve) 



Sn, compared with the transition density 



single-particle orbits below the N = 82 gap (ejjF < —7.77 
MeV) by a factor of twenty. Clearly the long tail associ- 
ated with the [3p 3 / 2 ] 2 configuration is a consequence of 
the weak binding of the 3p 3 / 2 quasineutron state. 

We speculate that many two-quasiparticle configura- 
tions including [3p 3 / 2 ] 2 and [2/ 7 / 2 ] 2 contribute to pro- 
duce the pair vibrational mode in 134 Sn as the transition 
density is larger by several times than those of the indi- 
vidual unperturbed two-quasiparticle excitations. Com- 
paring the transition density of [3p 3 / 2 ] 2 with that of the 
pair vibrational mode, we deduce that the long tail in 
the pair vibrational mode may be inherited from that of 
[3p 3 / 2 ] 2 . Two-quasineutron configurations in the contin- 
uum, including the next lying 3pi/ 2 unbound resonance 
state (located at positive energies around 0.2 MeV) may 
contribute also. Considering the orthogonality of the 
wave function of the pair vibrational mode to that of the 
pair rotational mode whose main component is [2p 7 / 2 ] 2 , 
we deduce that the [3p 3 / 2 ] 2 and [3pi/ 2 ] 2 configurations 



are the largest components in the pair vibrational mode. 

In Fig. EJa), we plot the pair-removal strength func- 
tion arising from unperturbed neutron two-quasiparticle 
excitations in 142 Sn. The low-lying peak seen at E = 3.3 
MeV is the [2/ 7 / 2 ] 2 two-quasiparticle configurations, cor- 
responding to creating two holes (i.e., removing two neu- 
trons) in the 2/ 7 / 2 orbit. There are also peaks associ- 
ated with the unperturbed [3p 3 / 2 ] 2 and [3p!/ 2 ] 2 config- 
urations, which are however hardly visible in this scale. 
It is seen that the strength of the pair vibrational mode 
in the full QRPA calculation is ten times larger than 
the unperturbed strength, indicating a collective effect. 
The transition densities of the two-quasiparticle configu- 
rations [2/ 7 / 2 ] 2 and [3p 3 / 2 ] 2 are shown in Fig.[6fb). Com- 
parison of the full and the unperturbed transition den- 
sities suggests that [2/ 7 / 2 ] 2 is one of the main configu- 
rations, but a significant mixing of configurations other 
than [2/ 7 / 2 ] 2 is also present. Since it is the pair-removal 
mode the configuration mixing involving those below the 



N = 82 shell gap can contribute, but the weakly bound 
and continuum orbits such as 3p 3 / 2 and ?>Pi/2 contribute 
very little to the pair removal mode as their occupation 
are small. Consequently the spatial extension of the tran- 
sition density of the pair-removal vibrational mode in 
142 Sn is smaller than that of the pair-addition vibrational 
mode in 134 Sn. 



D. Systematics 

It is interesting to see systematical behavior of the low- 
lying pair vibrational modes along the isotopic chain of 
Sn. Here we treat separately the pair-addition and the 
pair-removal modes of the pairing vibration. They are 
identified as the QRPA eigenmode having the largest 
pair-addition [-removal] strength among the peaks below 
E = 5 MeV. 

In Fig. [7] and Table U shown are the excitation energy 
and the pair-addition strength B(PadO) of the low-lying 
pair-addition vibrational modes. (We plot here the re- 
sults obtained with the DDDI-bare', the mix and the 
volume pairing interactions, but we shall discuss depen- 
dence on the pair interaction later.) We plot in Fig.[7Ja) 
also the calculated threshold energies S\ n and S^n of 
one- and two-neutron separations. A noticeable feature 
is that in the isotopes with A > 132 the low-lying pair 
vibrational modes are located near the neutron separa- 
tion energies. This is due to the sudden decrease of the 
separation energies in isotopes beyond the magic number 
N = 82. The minimum at A = 140 (N = 90) reflects the 
sub-shell closure of the 2/7/2 orbit. In Fig. [7] (b), large 
pair-addition strength _B(PadO) is seen in five isotopes 
with A = 132 - 140 (N = 82 - 90). 

Figure [8] shows the pair-addition transition density 
PiL=o ( r ) °f the low- lying pair-addition vibrational modes 
in the isotopes from A = 120 (a stable isotope) to 
A = 150. We find that the pair vibrational modes in 
the A = 132 — 140 isotopes have quite similar radial pro- 
files of the transition density i^^o(r)) and they all share 
the common character concerning the spatial extension 
reaching r ~ 15 fm. These features all indicate that the 
anomalous pair vibrational mode appears systematically 
in the region beyond the N = 82 magic number and up 
to N — 90. It is also evident that the anomalous pair vi- 
brational modes in A — 132 — 140 are distinct from those 
in the A — 120 — 130 isotopes closer to the stability line. 

Systematics of the pair-removal vibrational mode is 
shown in Fig. [S] and also in Table |TJ The excitation en- 
ergy and the pair-removal strength i?(Prm0) are plotted 
in Fig. |U[a) and (b), respectively. In the isotopes with 
A = 134 — 140, the low- lying pair-removal strength is 
negligible. However, the pair-removal vibrational mode 
having large strength emerges in the isotopes beyond 
A = 140 (iV = 90). It is found that the pair vibrational 
mode in these isotopes have the same character which 
we discussed in connection with the results for 142 Sn (cf. 
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Fig.Hd)). 
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3.138 


3.55 


1.180 


3.55 


0.003 


104 


5.101 


2.95 


1.023 


2.95 


0.037 


106 


6.498 


2.62 


0.595 


2.62 


0.183 


108 


7.857 


4.22 


0.225 


2.82 


0.370 


110 


8.932 


3.69 


0.245 


3.28 


0.405 


112 


9.872 


3.28 


0.198 


3.81 


0.435 


114 


9.798 


4.30 


0.820 


4.30 


0.448 


116 


9.333 


3.61 


0.800 


2.78 


0.488 


118 


8.855 


2.84 


0.493 


3.25 


0.888 


120 


8.274 


2.40 


0.402 


3.80 


0.984 


122 


7.851 


2.40 


0.136 


4.63 


0.920 


124 


7.374 


2.75 


0.031 


2.75 


0.699 


126 


6.462 


3.20 


0.010 


3.20 


0.625 


128 


4.991 


3.65 


0.004 


3.65 


0.557 


130 


2.851 


4.78 


0.005 


4.78 


0.498 


132 


3.202(r) 3.973(a) 


4.28 


2.747 


4.43 


0.458 


134 


3.607 


3.73 


3.160 


3.73 


0.000 


136 


5.662 


3.10 


3.683 


3.10 


0.007 


138 


6.048 


2.33 


4.615 


2.33 


0.058 


140 


5.157 


1.23 


2.929 


1.23 


2.376 


142 


11.079 


1.90 


0.256 


2.35 


2.800 


144 


16.088 


1.75 


0.182 


3.10 


1.983 


146 


19.869 


2.20 


0.119 


3.65 


1.563 


148 


22.639 


2.10 


0.082 


4.08 


1.338 


150 


24.481 


2.03 


0.071 


4.45 


1.077 



TABLE I: The pair-addition and -removal strengths 
B(Pad/rm0) of the ground-state transition, the excitation en- 
ergy E and the pair-addition strength B(Pad0) of the pair- 
addition vibrational mode, and the excitation energy E and 
the pair-removal strength _B(Pad0) of the pair-removal vibra- 
tional mode in Sn isotopes. The unit of the excitation energy 
is MeV. The ground-state transitions in the closed shell nu- 
clei 100 Sn and 132 Sn are evaluated in the QRPA, separately 
for the pair-addition and pair-removal transitions. 



E. Comparison with pair rotational ground-state 
transition 

The anomalous feature of the pair vibrational modes 
in A — 132 — 140 is clarified further by comparing with 
the pair transfer populating the ground state, i.e., the 
pairing rotation discussed in Section IHII 

Figure [TUta) show the ratio r = i?(Pad0;gs — > 
pv)/i?(Pad0; gs — > gs) of the pair-addition strength as- 
sociated with the low-lying pair vibrational mode and 
that with the pair rotational ground-state transfer. The 
strength ratio amounts to 60-90 % in the A = 132 — 140 
isotopes. It is compared with isotopes in the middle of 
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FIG. 7: (a)The excitation energy of the pair-addition vibrational mode. The one- and two-neutron separation energies are 
plotted with the dotted and the dashed lines, respectively, (b) The pair-addition strength B(Pad0) of the pair-addition 
vibrational mode. The diamonds connected with the solid line are the results obtained with the pairing interaction DDDI-bare' 
while the circles and the crosses are those with the mix and the volume pairing interactions, respectively. 




the N = 50 — 82 shell, where the ratio is far below 10 % 
(except in a few isotopes at the begging of the shell) . The 
small ratio in the isotopes A ~ 110 — 130 is in qualita- 
tive agreement with the well established observation that 
two-neutron transfer cross sections populating low-lying 
+ states is generally weak, for instance, less than 10% 
in stable Sn isotopes || [l2j. A known exception is the 
case where an excited + state having the character of 
the shape coexistence emerges as a consequence of sudden 
shape changes with the neutron number pj.l43|. The large 
strength ratio of 60-90 % in the A = 132— 140 isotopes is 
comparable with the cases of the shape transition/ shape 
coexistence. 

In Fig. UOf b) , we compare the pair-addition transition 

densities Pj^= (r) for the ground-state transitions and 
for the pair- vibrational transitions populating excited + 
states. The maximal values of the amplitudes are com- 
parable, but the transition density of the anomalous pair 
vibration mode extends more than that associated with 
the pair rotational ground-state transition. It reflects 
different microscopic structures of the pairing rotation 



and of the anomalous pairing vibration. As we discussed 
in connection with Fig. [5J the largest component of the 
pairing rotation may be [2/7/2] 2 , while the counter part 
of the anomalous pairing vibration may be [3p 3 / 2 ] 2 and 
[3p 2 /i] 2 - The 3p orbits have much smaller binding energy 
and thus longer tail in the wave functions than those of 
2/7/2- The difference in the spatial extension of the tran- 
sition densities may be explained in this context. 



The ratio r' — £>(Prm0;gs — > pv)/B(Prm0; gs — > gs) 
of the strengths of the pair-removal vibrational mode and 
of the pairing rotation is shown in Fig. HOT c). The ratio 
in 14 °- 144 Sn is 30 — 40%, but it is relatively small if com- 
pared to the ratio 60 — 90% of the pair addition strengths 
in 132_140 Sn, and it decreases to small values in A > 146. 
Note that in the isotopes beyond A = 140 the strengths 
of the pair rotational ground-state transfer is more dras- 
tically increased than the strengths of the pair-removal 
vibrational mode. 
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FIG. 9: (a) The excitation energy and (b) the pair-removal strength B(PrmO) of the pair-removal vibrational mode. See also 
the caption of Fig. [7) 




FIG. 10: (a) The ratio of the pair-addition strength B(PadO; gs — Y pv) of the low-lying pair-addition vibrational mode and the 
strength B(PadO;gs — > gs) of the ground-state transfer, (b) The pair-addition transition density r 2 P^=o( r ) associated with 
the anomalous pair vibrational mode ( solid curve) and that for the ground-state transfer ( dashed curve) in 134 Sn. (c) The 
ratio of the pair-removal strength B(PrmO; gs — > pv) of the low- lying pair-removal mode and the ground-state transfer strength 
B(PadO;gs -> gs). 



V. SENSITIVITY TO DENSITY-DEPENDENT 
PAIRING 

It is interesting to examine the sensitivity of the pai r 
transfers to different effective pairing interactions pll [23j . 
For this purpose we have performed the HFB+QRPA cal- 
culations with the three different versions of the DDDI 



described in Section III Al Note that the three interac- 
tions are designed to give approximately the same values 
of the average neutron pair gap A uv (cf. Fig.QJb)). Com- 
paring the three interactions, we can examine effects of 
different density-dependences of the DDDI while keeping 
approximately the same average pair gap. 

The influence of the density dependence on the pair- 
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FIG. 11: (a) The pair-addition strength function iSpado(-E) for neutrons in 134 Sn, calculated with the three effective pairing 
interactions; the DDDI-bare', the mix, and the volume types. The arrows indicate the peak heights, (b) The transition density 

r 2 P>^2 (r) of the low-lying pair-addition vibrational mode in 134 Sn, calculated with the three effective pairing interactions; the 
DDDI-bare', mix, and volume types. In the inset we show the transition density P}i^ (r) without the volume element. 



addition and -removal strengths i?(Pad/rmO) of the 
ground-state transition is shown in Fig. [IJa) . Although 
the influence is not significant for A = 100 — 132, we see 
some sensitivity (but not very strong one) to the density 
dependence for neutron-rich isotopes with A > 132. 

Concerning the pair vibrational modes, in contrast, the 
anomalous pairing vibration in A = 132 — 140 has sig- 
nificant sensitivity to the density-dependence of the ef- 
fective pairing interactions. As an example, we show in 
Fig. IrTTa') the pair-addition strength function Sp a do(E) 
in 134 Sn calculated with the three different DDDI's. It is 
noticed that the basic structure of the strength function 
is the same, and the anomalous pair vibrational mode 
commonly exists at around E ps 3.5 MeV. More impor- 
tantly it is seen that the height of the pair vibrational 
peak depends strongly on the different pairing interac- 
tion by about a factor of 1.5. The difference is larger 
than that found in the ground-state transition. The same 
effect is seen commonly at A — 132 — 138 as shown in 
Fig-EJb). On the other hand, the influence of the density- 
dependence on the pair-addition vibrational in the iso- 
topes A < 132 mode is small. 

We can understand the influence of the density- 
dependence on the anomalous pair vibrational mode in 

terms of the transition density -P/2=o( r )' wnose depen- 
dence on the pair interactions is shown in Fig. Illf b). 

It is seen in this figure that the amplitude Pil= (r) m 
the external region r <; i? rms outside the surface shows 
the same trend as in the pair- addition strength S(Pad0); 
the DDDI-bare' > mix > volume pairing interactions. 
We can then relate this trend to the features of density- 
dependent interaction strength V n (r), which is, in the 
external region, the largest for the DDDI-bare' and the 
smallest for the volume pairing, reflecting the density- 
dependent interaction strength at low densities, (vg = 
-458, -292, -195 MeV fm 3 and see also Fig.l of Ref.[23| 
). Since the anomalous pair vibrational mode at A — 
132 — 140 has the transition density extending far out- 



side the surface, the large difference in the interaction 
strength V n (r) in the exterior is effective to this mode. 
However, for the usual pair vibrational modes whose 
transition density does not reach very far, the sensitivity 
becomes small. 

Finally we remark on related works [U [23[ discussing 
the sensitivity to the density-dependence of the effec- 
tive pairing interaction. In Ref. [2l|, the pair-transfer 
strengths and the transition density of the pair vibra- 
tional modes are analyzed for 124 Sn and 136 Sn by using 
the Skyrme-HFB + QRPA model, which is similar to the 
model adopted here. In the case of 136 Sn, the low-lying 
pair vibrational mode around E ~ 3 MeV emerges in 
two of their calculations (the rj = 0.65 and 0.35 cases of 
Ref. plj ) . in qualitative agreement with our results. The 

shape of the pair-addition transition density P^^} Q {r) is 
also similar as seen from the comparison of the inset of 
Fig. [Tl7b) of the present paper and Fig. 8 of Ref. [21]] . 
On the other hand, in another calculation adopting the 
DDDI with the strongest density dependence (the r\ = 1 
case in Ref. (2l| ) . the low- lying pair vibrational mode 
is fragmented into three peaks (cf. Fig. 7 of Ref . [2ll] ) . 
suggesting a complex sensitivity to the pairing interac- 
tions. We note that the interaction strength vo = —670 
MeVfm 3 chosen in the rj = 1 case is too large since it 
leads to the positive scattering length, implying unphys- 
ical existence of a bound state for the 5 = neutron 
pair. In the present investigation, on the contrary, we 
constrain the upper limit of \vq\ as \vo\ < 458 MeV fm 3 
by the scattering length a = —18.5 MeV of the 1 S I chan- 
nel of the bare nuclear force. For our pairing interactions 
determined under this constraint, the anomalous pair vi- 
brational mode stably emerges as a single-peak around 
E ~ 3 - 4 MeV in 132 " 140 Sn. 

It is discussed in Ref. [23] that the quadrupole pair- 
addition transfer populating the first 2 + states in 
134 Sn and heavier isotopes is sensitive to the density- 
dependence of the pairing interaction. In that case the 
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ratio of the pair-addition strengths between the volume 
pairing and the DDDI-bare' amounts to approximately 
a factor of two, which is slightly larger than the sensi- 
tivity (a factor of ~ 1.5) of the pair vibrational mode in 
132-I40g n . The origin of the different sensitivity is not 
clear at present, but it may be related to the fact that 
the monopolc pair transition density has larger ampli- 
tude also in the internal region of the nucleus than that 
of the quadrupole mode. 

VI. CONCLUSIONS 

We have investigated microscopically monopole two- 
neutron transfer modes in heavy-mass supcrfluid nuclei 
using the Skyrme-Hartree-Fock-Bogoliubov mean-field 
model and the continuum quasiparticle random phase 
approximation. Emphases have been put on the pair vi- 
brational modes populating low-lying collective + states 
and on the pair rotational transitions connecting the 
ground states of the neighboring AN = ±2 isotopes. 
We have performed systematic numerical analysis for the 
even-even Sn isotopes ranging from A = 100 to A = 150, 
and found the following new features of the two-neutron 
transfer modes emerging in neutron-rich Sn isotopes be- 
yond the N = 82 shell closure (A > 132). 

The calculation predicts a novel type of pair vibra- 
tional mode in the isotopes 132 ~ 140 Sn. It is a monopole 
vibrational mode characterized by intense transition from 
the ground state to a low-lying + state in the neighbor- 
ing N + 2 isotope via the two-neutron addition transfer. 
The corresponding + states are predicted to emerge near 
the threshold energy of the one-neutron separation, and 
to form a narrow resonance in case when it is located 
above the separation energy. An anomalous feature of 



this pair vibrational mode and a marked difference from 
the pairing vibration in nuclei near the stability line is 
that the pair transition density exhibits a large amplitude 
in the region outer than the surface, and also a long tail 
extending up to r <~ 15 fm. We expect that the differ- 
ence may manifest itself in reaction observables such as 
the cross section of the (t,p) reaction although quantita- 
tive analysis remains for future investigations. As a mi- 
croscopic origin of the anomalous pair vibrational mode, 
the weakly bound or resonant 3p 3 / 2 and 3pi/2 orbits are 
suggested to play key roles. 

In very neutron-rich isotopes with A > 140 (N > 90), 
enhancement of the two- neutron transfer strength is pre- 
dicted for the transitions between the ground states. 
This is a natural consequence of the spatially extended 
pair field originating from the weak binding of neutrons. 
Particularly, occupation of the weakly bound 3p 3/ / 2 and 
3pi/2 orbits gives the sudden increase in the strength for 
A > 140. In this perspective, the anomalous pair vi- 
brational mode in 132 ~ 140 Sn can be linked to the strong 
ground-state transfer in A > 140, and it can be re- 
garded as a precursor of the enhanced pair transfer which 
emerges in weakly bound neutron-rich nuclei close to the 
drip-line. 
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